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ABSTRACT 

In this paper we present multi-band period-luminosity (P-L) relations for fundamental mode 
Cepheids in the SMC. The optical U/-band mean magnitudes for these SMC Cepheids were taken 
from the third phase of the Optical Gravitational Leasing Experiment (OGLE-III) catalog. We also 
matched the OGLE-III SMC Cepheids to 2MASS and SAGE-SMC catalog to derive mean magnitudes 
in the JiLAT-bands and the four Spitzer IRAC bands, respectively. All photometry was corrected for 
extinction by adopting the Zaritsky’s extinction map. Cepheids with periods smaller than ~ 2.5 days 
were removed from the sample. In addition to the extinction corrected P-L relations in nine filters 
from optical to infrared, we also derived the extinction-free Wesenheit function for these Cepheids. We 
tested the nonlinearity of these SMC P-L relations (except the 8.0/im-band P-L relation) at 10 days: 
none of the P-L relations show statistically significant evidence of nonlinearity. When compared to the 
P-L relations in the LMC, the t-test results revealed that there is a difference between the SMC/LMC 
P-L slopes only in the V- and J-band. Further, we found excellent agreement between the SMC/LMC 
Wesenheit P-L slope. The difference in LMC and SMC Period-Wesenheit relation LMC and SMC zero 
points was found to be A/r = 0.483 ± 0.015 mag. This amounts to a difference in distance modulus 
between the LMC and SMC. 

Subject headings: Magellanic Clouds — stars: variables: Cepheids — distance scale 


1. INTRODUCTION 

The period-luminosity relation (the Leavitt Law, here¬ 
after P-L relation) for Small Magellanic Cloud (SMC) 
Cepheids was first presented in Leavitt & Pickering 
(1912). Since then, the SMC P-L relations have been 
derived from optical to near-infrared (for examples, see 
Arp I960; Payne-Gaposchkin 1965; Wayman 1984; Welch 
& Madore 1984; Visvanathan 1985; Welch et al. 1985; 
Caldwell & Coulson 1986; Laney & Stobie 1986; Math- 
ewson et al. 1986; Welch et al. 1987; Caldwell & Laney 
1991; Smith et al. 1992; Laney & Stobie 1994; Nemec 
1994; Sharpee et al. 2002, and references therein) based 
on relatively small number of SMC Cepheids. In 1999, 
the second phase of the Optical Gravitational Leasing 
Experiment (hereafter OGLE-H) released a catalog that 
contained more than 450 fundamental mode Cepheids 
located in ^ 2.4 square degree at the center of SMC 
(Udalski et al. 1999a). A number of SMC P-L relations 
have been derived in literature based on these OGLE- 
11 SMC Cepheids (Udalski et al. 1999b; Groenewegen 
2000; Storm et al. 2004; Tammann et al. 2008; Sandage 
et al. 2009; Bono et al. 2010). Independently, the EROS 
(Experience de Recherche d’Objets Sombres) Collabora¬ 
tion also derived the P-L relations based on a large num¬ 
ber of SMC Cepheids in customized filters (Sasselov et 
al. 1997; Bauer et al. 1999; Marquette 1999). 

In 2010, a catalog for an even larger number of 
Cepheids in SMC was released from the the third phase 
of OGLE operation (hereafter OGLE-Hl, see Soszyhski 
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et al. 2010). Compared to OGLE-H, more than 2600 
fundamental mode SMC Cepheids were included in the 
OGLE-Hl SMC catalog as a result of larger survey area 
(Soszyhski et al. 2010). Although the U/-band SMC P- 
L relations were also derived in Soszyhski et al. (2010), 
these P-L relations were not corrected for extinction. 
Subsequently, Subramanianfe Subramaniam (2015) used 
the OGLE-Hl UJ-band photometry for Cepheids to in¬ 
vestigate the spatial structure of SMC. In terms of the 
near infrared JJJJJ-band P-L relations, Matsunaga et al. 
(2011) presented preliminary SMC P-L relations in three 
period bins by matching the OGLE-Hl SMC Cepheids to 
the single-epoch Magellanic Clouds point source catalogs 
(Kato et al. 2007) based on the Infrared Survey Facility 
(IRSF) observations. Inno et al. (2013) further combined 
these IRSF measurements with OGLE-III UJ-band pho¬ 
tometry to derive the period-Wesenheit relations in vari¬ 
ous combinations. It is expected that more near infrared 
data for the SMC Cepheids will be available from the 
VISTA survey of the Magellanic Clouds System Project 
(VMC Cioni et al. 2011) in the near future. 

Using the OGLE-III LMC Cepheid catalog (Soszyhski 
et al. 2008), Ngeow et al. (2009, hereafter Paper I) de¬ 
rived the extinction corrected P-L relations in VIJHK 
and the four Spitzer IRAC bands. In this work, we ex¬ 
tend our investigation and derive the extinction corrected 
multi-band P-L relations based on the OGLE-III SMC 
Cepheids using the catalog from Soszyhski et al. (2010). 
It may be noted that the SMC IRAC bands P-L relations 
were derived in Ngeow & Kanbur (2010) based on single 
epoch Spitzer data. In this work the IRAC bands P- 
L relations are updated using the available photometry 
up to three epochs. The SMC P-L relation is particu¬ 
larly important in distance scale and stellar pulsation 
work, such as constraining the theoretical predictions 
(see Bono et al. 2010). This is because the metallicity 
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of SMC is 12 + \og{0/H) = 7.98 dex, which is similar 
or comparable to other local dwarf galaxies (such as IC 
1613, 7.86 dex; WLM, 7.74 dex; Sextans A, 7.49 dex; 
Sextans B, 7.56 dex; Pegasus, 7.92 dex; Leo A, 7.38 dex; 
see Sakai et al. 2004; Tammann et al. 2011). 

2. DATA AND EXTINCTION CORRECTION 

Mean C/-band magnitudes and periods for 2626 funda¬ 
mental mode SMC Cepheids were taken from Soszyhski 
et al. (2010). The Wesenheit function, W = I— 1.55(y— 
/), was also calculated from the mean C/-band mag¬ 
nitudes (if the C-band mean magnitude is available). 
The OGLE-111 SMC Cepheids were also matched to the 
2MASS point source catalog (Cutri et al. 2003; Skrutskie 
et al. 2006), using a search radius of 2". Mean separa¬ 
tion of the 2281 matched 2MASS sources is 0.225", with 
a dispersion of 0.264" (only 60 matched 2MASS sources 
have separation greater than 1"). Random-phase correc¬ 
tions as described in Soszyhski et al. (2005) were applied 
to 2MASS photometry to derive the mean JHK magni¬ 
tudes, using the scaling between /-band amplitudes and 
the Ji/AT-band amplitudes. Finally, up to three epochs 
of the IRAC band photometry, based on publicly released 
SAGE-SMC (Surveying the Agents of Galaxy Evolution 
in the Tidally Disrupted, Low-Metallicity Small Magel¬ 
lanic Cloud, Gordon & SAGE-SMG Spitzer Legacy Team 
2010) data, were downloaded from the Spitzer Science 
Genter. As in Paper 1 and Ngeow & Kanbur (2010), 
we adopted the SAGE-SMC archival data (version S14 
and later, delivered on 2010 September 30) in this work. 
A search radius of 2" was used to match the OGLE-111 
SMG Cepheids and the sources in SAGE-SMC archival 
data. The number of matched sources and the corre¬ 
sponding mean separations are summarized in Table 1 
for the SAGE-SMC Epoch 0, 1 and 2 data. Intensity 
means were calculated using the three epochs data (when 
available) for each matched Cepheids in the IRAC bands. 

As in Paper 1, extinction for each OGLE-111 SMG 
Gepheid was estimated using the Zaritsky et al. (2002) 
extinction map. For a given input location of SMG 
Gepheids, this extinction map returns the extinction 
in E-band {Ay), measured from the cool stars only. 
In case the extinction maps did not return any ex¬ 
tinction values for a given Gepheid, a mean value of 
Ay = 0.18 was adopted. Extinctions in other bands were 
scaled using the following total-to-selective extinction 
coefficient. Ry^ /, j, Je. S.S/xm, I.S/xm, S.S/xm, S.O^m — 
{3.24, 1.96, 0.95, 0.59, 0.39, 0.17, 0.12, 0.08, 0.05}. 
Following Paper 1, the total-to-selective extinction co¬ 
efficients in EJ bands are adopted from Udalski et al. 
(1999b), and in other bands these values are calcu¬ 
lated based on the extinction law from Gardelli et al. 
(1989). Our values are slightly different from the value 
of Ry = 3.1 adopted in Zaritsky et al. (2002) extinction 
map (Zaritsky 1999), and almost identical to those used 
by Fouque et al. (2007) in VIJHK bands. 

3. THE PERIOD-LUMINOSITY RELATIONS 

Extinction corrected mean magnitudes in each band 
were used to derive the corresponding P-L relations. The 
resulting P-L relations are presented in Figures 1 and 
2. We restricted our Gepheid sample to period range 
0.4 < log(P) < 1.515 only. Soszyhski et al. (2010) rec¬ 
ommended against using the 17 brightest Gepheids listed 


TABLE 1 

Summary of the Matched SAGE-SMC Archival Data. 


Band 

-^match 

< A >=■ 

Std. Dev.b 

Fraction within 1"^ 


Epoch 0 

3.6/^m 

1370 

0.314 

0.260 

97.23% 

4.5/^m 

1367 

0.299 

0.254 

97.59% 

5.8/im 

467 

0.287 

0.310 

95.50% 

8.0/^m 

361 

0.281 

0.317 

96.12% 

Epoch 1 

3.6/^m 

2580 

0.272 

0.241 

98.22% 

4.5/i.m 

2557 

0.271 

0.242 

98.20% 

5.8/i.m 

702 

0.276 

0.299 

96.30% 

8.0/^m 

406 

0.295 

0.313 

95.57% 

Epoch 2 

3.6/Tm 

2557 

0.309 

0.242 

97.89% 

4.5/im 

2545 

0.308 

0.238 

97.92% 

5.8/im 

686 

0.318 

0.277 

96.65% 

S.O/rm 

366 

0.334 

0.298 

95.90% 


^ A is the separation, in arcsecond, between the matched SAGE-SMC 
Archival sources and the OGLE-III SMC Cepheids. 

^ The standard deviation of the mean. 

^ Fraction of matched SAGE-SMC Archival sources within l" radius 
from the OGLE-III SMC Cepheids. 

in the /-band catalog for absolute calibration of bright¬ 
ness. The reason is that these Cepheids were not pro¬ 
cessed with the standard OGLE Difference Image Anal¬ 
ysis (DIA) pipeline as these stars saturate in the refer¬ 
ence frames. The shortest period of these 17 Cepheids 
is log(P) ~ 1.515, therefore we adopted an upper limit 
of the period at log(P) = 1.515. Furthermore, Figures 
1 and 2 suggest that the slopes of the P-L relations for 
Cepheids with log(P) > 1.9 become flatter in several 
bands, as they belong to a sub-class of Cepheids - the 
ultra-long period Cepheids (Bird et al. 2009). The SMC 
short period Cepheids are known to exhibit a change 
in the slope of the P-L relation at a specific period. 
This break period at log(P) ~ 0.4 has been seen in 
optical data (Bauer et al. 1999; Udalski et al. 1999b; 
Sharpee et al. 2002; Sandage et al. 2009; Soszyhski et 
al. 2010; Bhardwaj et al. 2014) and extended to mid- 
infrared (Ngeow & Kanbur 2010). However, Tammann 
et al. (2011) suggested the break period might occur 
at log(P) = 0.55, while Subramanian & Subramaniam 
(2015) found that it could happen at log(P) = 0.47. Nev¬ 
ertheless, we adopted log(P) = 0.4 as the short period 
end of our sample. Since our goal is to derive the P-L 
relations for distance scale applications using the long 
period Cepheids, detailed investigation into these SMC 
Cepheids with log(P) < 0.4 will be presented elsewhere. 
Finally, we removed two additional Cepheids (OGLE- 
SMC-CEP-1157 and OGLE-SMC-CEP-3212), because 
they were marked in the remarks. txt file as their pho¬ 
tometry was derived from the DoPHOT package rather 
than the DIA method. 

Inspecting Figures 1 and 2 suggests that a period cut 
is needed for K-, 5.8 and 8.0/im-band P-L relations, after 
removing log(P) < 0.4 Cepheids, to avoid the influence of 
incompleteness bias at the short-period end. This incom¬ 
pleteness bias is due to the well-known Malmquist bias 
as discussed, for example, in Sandage (1988). We per¬ 
formed the following steps to determine the appropriate 
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Fig. 1.— The extinction corrected SMC P-L relations in VI bands (left panels) and JHK bands (right panels). The extinction-free 
Wesenheit function is also included in lower-left panel. The (green) filled squares represent the excluded Cepheids at which the OGLE- 
III photometry was based on the DoPHQT package (see text for more details). The dashed lines in right panels represent the 2MASS Scr 
sensitivity adopted from Cutri et al. (2003). 
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Fig. 2. — The extinction corrected SMC P-L relations in 
IRAC bands for all the matched sources. The (green) filled 
squares represent the excluded Cepheids at which the OGLE- 
III photometry was based on the DoPHOT package (see text for 
more details). The horizontal dashed lines represent the faint 
limits in each bands, taken from SAGE-SMC document (See 
http://data.spitzer.caltech.edu/popular/sage-sinc/20100930 
.enhanced/document s/sage--sinc_delivery_seplO .pdf). 

period cuts in the K-, 5.8 and 8.0/rm bands, respectively. 

1. We first adopted an initial period cut at log(Pcut) = 
0.4 and removed Cepheids having periods smaller 
than this value. 

2. We fitted a P-L relation to the remaining sample of 
Cepheids with an iterative 2.5(T clipping algorithm 
to derive the P-L slopes (and its associated error) 
for this sample of Cepheids. 

3. We repeated step 1 and 2 with a binsize in 


AlogP = 0.01, up to a maximum value of 
log(Pcut) = 1-01. We picked this binsize such that 
the parameter space of log(Pcut) can be properly 
sampled. Upper panels of Figure 3 display the fit¬ 
ted P-L slopes as a function of adopted period cuts 
in K, 5.8 and 8.0^m bands. 

4. We then plotted the distributions of the P-L slopes 
given in upper panels of Figure 3 and present these 
distributions as histograms in the lower panels. 

5. Based on these histograms, we determined the 
mode of the histogram as an indicator that the 
P-L slopes begin to stabilize (i.e. without the in¬ 
fluence of incompleteness at the short period-end) 
at a given period cut. The values of the mode are 
shown as horizontal and vertical dashed (red) line 
in upper and lower panels of Figure 3, respectively. 

6. We calculated the absolute deviation between the 
P-L slopes given in the upper panels of Figure 3 and 
the modal value determined from the histograms 
(i.e. the absolute deviation between the horizon¬ 
tal dashed line and individual opened circles in the 
upper panels of Figure 3). 

7. The P-L slope that returned the smallest value of 
the absolute deviation, based on step 6, is adopted 
as the final log(Pcut) as indicated by a downward 
vertical arrow in the upper panels of Figure 3. In 
case more than one P-L slope resulted in the small¬ 
est value of absolute deviation, we adopted the one 
with minimum value of log(Pcut)- 

Hence, the final adopted log(Pcut) in K-, 5.8/rm- and 
8.0^m-band are 0.51, 0.68 and 0.91, respectively. 

Following the OGLE team (Udalski et al. 1999b; 
Soszyhski et al. 2010), obvious outliers of the P-L re¬ 
lations, as displayed in Figures 1 and 2, are removed 
using an iterative 2.5cr clipping algorithm (cf. Paper I), 
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Fig. 3. — The top panels displayed slopes of the fitted P-L relations for data points with log{P) > log(Pcut) as a function of the adopted 
period cut. Histograms of the distribution of these fitted P-L slopes were given in the bottom panels. The (red) dashed-lines indicate the 
mode of the fitted P-L slopes based on the histograms. For 8.0/rm-band P-L slopes, we only derive the mode value from slopes that are 
smaller than —3.15. The vertical arrows represent the adopted log(Pcut) in a given bands. Errors on the P-L slopes are standard errors 
calculated from the ordinary least squares (OSL) regression method. 


TABLE 2 

Multi-Band P-L Relations for SMC Cepheids. 


Band 

P-L Slope 

P-L ZP 

P-L Dispersion (( 7 ) 

V... 

-2.660 ± 0.040 

17.606 ± 0.028 

0.275 

I... 

-2.918 ± 0.031 

17.127 ±0.022 

0.214 

J... 

-3.052 ± 0.026 

16.773 ± 0.019 

0.182 

H... 

-3.157 ±0.025 

16.530 ±0.018 

0.174 

K... 

-3.213 ±0.032 

16.514 ±0.025 

0.174 

3.6//m 

-3.220 ±0.021 

16.433 ±0.015 

0.146 

4.5/i.m 

-3.184 ±0.022 

16.375 ±0.016 

0.155 

5.8/i.m 

-3.201 ±0.039 

16.362 ± 0.036 

0.143 

8.0//m 

-3.268 ± 0.087 

16.426 ±0.097 

0.156 

w... 

-3.314 ±0.020 

16.375 ±0.014 

0.137 


where a represents the dispersion of the P-L relation in 
each iteration. These outliers could be caused by a va¬ 
riety of reasons including, but not limited to, blending 
of nearby sources, mis-matching of the sources and mis- 
identification as classical Cepheids in the OGLE-III cata¬ 
log. For example, Ngeow & Kanbur (2010) discussed the 
outliers found in the IRAC bands. Since the goal of this 
work is to derive the P-L relations, we do not investigate 
the reasons or nature of these outliers in detail. Never¬ 
theless, it is clear that these outliers should be removed in 
fitting the P-L relations. Figure 4 shows the correlation 
of residuals of P-L relations in two bands. This demon¬ 
strates that the majority of these outliers are present in 
either or both bands. We have also tried to fit the P-L 
relations, including the outliers, using the robust regres¬ 
sion technique (a regression technique that is robust to 
outliers presented in the sample. Street et al. 1988; Du¬ 
Mouchel & O’Brien 1989). The differences obtained in 
the P-L slopes and zero points using robust regression 
and our iterative 2.5(t clipping algorithm do not exceed 
0.015 and 0.012, respectively, in any given band. 

Our final derived P-L relations in various bands for 
the SMC Cepheids are summarized in Table 2. Slopes of 
the four IRAC bands and the P/VF-band P-L relations 
given in Table 2 are in good agreement with previous 
determinations as reported in Ngeow & Kanbur (2010), 
Soszyhski et al. (2010) and Subramanian & Subrama- 
niam (2015). 

3.1. Comparison with P-L Relations from OGLE-II 



-2 -1 0 1 2 
P-L residuals in band j^l 


Fig. 4.— Correlations of residuals of P-L relations in two bands. 
Crosses are the Cepheids remained in the sample after applying 
the iterative 2.5(7 clipping algorithm, while open (red) squares rep¬ 
resent the rejected outliers in the P-L relations (in either bands). 
The dashed boxes represent the 2.5(7 boundaries to reject the out¬ 
liers, where <j is the dispersion of the P-L relation in a given band. 

The multi-band P-L relations based on the OGLE-111 
SMG Gepheids are compared to the P-L relations de¬ 
rived from the OGLE-11 catalog in Table 3. As in Pa¬ 
per 1, we applied the t-test to test the consistency of 
P-L slopes derived here to the published values. The 
calculated T-values and the corresponding p-values (the 
probability), evaluated based on the t-distribution with 
a = 0.05 (where a is the adopted significance level), are 
listed in the last two columns of Table 3. The expected 
t-value from the t-distribution, which only depends on 
a and the degrees of freedom (z/), is 1.968 for v = 300 
or 1.962 for v = 1400 (which roughly covers the range 
of V in our test). The null hypothesis is that the P-L 
slopes given in Table 2 are the same as the P-L slopes 
based on OGLE-II Gepheids. This can be rejected if 
|T| > t, where t ~ 1.96 (or equivalently, p-value smaller 
than 0.05). From Table 3, the P-L slopes in all bands 
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TABLE 3 

Comparison of SMC P-L Relations. 


Ref. 

P-L Slope 

P-L ZP 

N 

m 

p-value 

1 

-2.660 ± 0.040 

y-band 
17.606 ± 0.028 

912 



2 

-2.572 ± 0.042 

17.480 ± 0.032 

466 

1.495 

0.135 

3 

-2.573 ±0.041 

17.492 ± 0.032 

464 

1.492 

0.136 

4 

-2.588 ±0.045 

17.530 ± 0.035=' 

460 

1.191 

0.234 

5 

-2.590 ±0.047 

17.600 ±0.019 

488 

1.145 

0.253 

1 

-2.918 ±0.031 

/-band 

17.127 it 0.022 

918 



2 

-2.857 ±0.033 

17.039 ± 0.025 

488 

1.336 

0.182 

3 

-2.843 ± 0.033 

17.052 ± 0.025 

487 

1.638 

0.102 

4 

-2.862 ±0.035 

17.083 ± 0.027= 

462 

1.192 

0.233 

5 

-2.865 ±0.036 

17.117 ±0.014 

488 

1.122 

0.262 

1 

-3.052 ±0.026 

J-band 

16.773 it 0.019 

883 



6 

-3.037 ±0.034 

16.771 ± 0.027 

418 

0.359 

0.720 

1 

-3.157 ±0.025 

/1-band 
16.530 it 0.018 

875 



6 

-3.160 ±0.032 

16.475 ± 0.025 

414 

0.075 

0.940 

1 

-3.213 ±0.032 

/^-band 
16.514 it 0.025 

627 



6 

-3.212 ±0.033 

16.494 ± 0.026 

418 

0.022 

0.983 

1 

-3.314 ±0.020 

IT-band 
16.375 ± 0.014 

909 



2 

-3.303 ± 0.022 

16.345 ± 0.017 

469 

0.369 

0.712 

3 

-3.310 ±0.020 

16.387 ±0.016 

463 

0.139 

0.890 

7 

-3.300 ±0.021 

16.381 ± 0.016 

446 

0.473 

0.636 


Note. — Reference: (1). this work; (2) Udalski et al. (1999b); (3). 
same as (2) but updated in ftp://sirius.astrouw.edu.pl/ogle/ 
ogle2/var_stars/smc/cep/catalog/README.PL; (4) Sandage et al. (2009); 
(5) Storm et al. (2004); (6) Groenewegen (2000); (7) Groenewegen 
(2000), with 2.5o'-clipping. 

^ By adopting fisMC — 18.93 mag. 

are consistent with those derived from OGLE-II catalogs. 
We did not consider the updated SMC P-L slopes given 
in Udalski (2000) as the author assumed the P-L slopes 
are the same in laoth LMC and SMC. For U/-band P-L 
relations, we note that different photometric reduction 
packages were used in OGLE-II (DoPHOT, see Udalski et 
al. 1998; Udalski 2000) and OGLE-III (DIA, see Udalski 
et al. 2008). For common Gepheids with log(P) > 0.4 
in OGLE-II and OGLE-III catalogs, the averaged differ¬ 
ence in mean magnitudes (i.e. OGLE-II minus OGLE- 
III) is ^ 0.009 mag and ~ 0.004 mag in V- and /-band, 
respectively. After including extinction corrections, the 
averaged differences increased to ^ —0.093 mag in the 
U-band and ~ —0.057 mag in the /-band. In addition, 
these differences do not show any dependency on the pul¬ 
sation period. 

4. TESTING EOR NONLINEAR P-L RELATIONS AT 10 
DAYS 

The LMC P-L relation is known to exhibit a break at 
a period around 10 days, as shown in Paper I (and refer¬ 
ences therein) and Bhardwaj et al. (2015) with rigorous 
statistical tests. Sandage et al. (2009) proposed that the 
same break should occur for the SMC P-L relations in 
optical bands based on an analysis of the period-color 
relations. However their derived SMC P-L relations do 
not exhibit a break at 10 days. In the Ji/AT-bands, Mat- 
sunaga et al. (2011) did not find any significant break 
at 10 days using the single-epoch data taken from the 
IRSF. In contrast, studies done in Bono et al. (2010) and 
Carcfa-Varela et al. (2013) suggested that the SMC P-L 
relations should be nonlinear. In this section we exam¬ 
ine the nonlinearity of the SMC P-L relations based on 


our data presented in previous section, where nonlinear¬ 
ity refers to an underlying P-L relation can be separated 
to two P-L relations with a break period at 10 days. Ta¬ 
ble 4 presents the fitted SMC P-L relations for long and 
short period Cepheids separated at 10 days. 

As pointed out by Ngeow & Kanbur (2006), statistical 
tests are needed to find out the nonlinearity of P-L rela¬ 
tion. We apply the F-test (as in Paper I) and a random 
walk method in this work to examine the nonlinearity 
of multi-band SMC P-L relations with a break period 
at 10 days. The F-test requires independent and identi¬ 
cally distributed (iid) variables, as well as homoscedas- 
ticity and normality of residuals. Since observations of 
one Cepheid are independent of observations of other 
Cepheids, the iid assumption is satisfied. We next ex¬ 
amine the homoscedasticity assumption (i.e. constant 
variance) of the P-L residuals. The left panel of Figure 
5 displays the residuals of P-L relation: these do not ex¬ 
hibit any obvious trends that violate the homoscedastic¬ 
ity assumption. Further, the average of the P-L residuals 
are found to be consistent with zero in all bands. Finally, 
we tested the assumption of normality of the P-L resid¬ 
uals. The middle and right panels of Figure 5 show the 
histogram and the quantile-quantile (qq) plot of the P-L 
residuals, respectively. The distribution of residuals can 
be represented with a Gaussian function, as indicated 
by a dashed curve in middle panel of Figure 5. The qq 
plot is a common diagnostic tool to evaluate the normal¬ 
ity assumption: the quantiles of the data should fall in 
a diagonal straight line (i.e. the y = x line). Our qq 
plot demonstrates that the majority of the P-L residuals 
follow a normal distribution (except few points at the 
extreme ends of qq plot). Even though in Figure 5 we 
only showed the results of P-L residuals from the We- 
senheit function, P-L residuals in other bands all look 
very similar to this Figure. Normality of the residuals is 
important in assuring that the distribution of the F-test 
statistic under the null hypothesis of linearity follows the 
F distribution. However we can also use a permutation 
method (as describe below, Kim et al. 2000) to gener¬ 
ate the distribution of the F statistic in a way that is 
independent of the distribution of the residuals. 

We fit a single regression line to the actual data and 
then fit two lines with a break point at a period of 
10 days, followed by calculating the F-test statistic, Fq. 
Details of the F-test can be found in Paper I and ref¬ 
erence therein, and will not be repeated here. In short, 
a Fq value that is greater than ~ 3 would indicate that 
the underlying P-L relation is nonlinear. If we have N 
data points we can calculate N residuals from the sin¬ 
gle straight line fit, (where i = Next we 

randomly permute the residuals = ej and add the per¬ 
muted residuals back to the fitted data from the single 
straight line fit. This generates one iteration of pseudo¬ 
data. Note the data point still has the original period. 
Now with this pseudo-data, we fit a straight line and then 
two lines with the original break point and calculate the 
F-test statistic, Fi. We repeat this process 10,000 times 
to obtain 10,000 independent F-test statistics. We then 
find the proportion of the Fi that are greater than the 
observed Fq . This proportion gives the p-value, or signif¬ 
icance, of the observed Fo-test statistic without assum¬ 
ing normality. Results of the F-test statistic, Fq and 
the p-value calculated using the permutation method are 



6 


Ngeow et al. 


TABLE 4 

SMC P-L Relations Separated at 10 Days. 


Band 

P-L Slopes 

P-L 

ZPs 


Ns 

P-L Slopes 

P-L 

ZPl 


Nl 

u... 

-2.634 

± 

0.061 

17.592 

± 

0.039 

0.269 

821 

-2.453 

± 

0.244 

17.345 

± 

0.290 

0.328 

91 

I... 

-2.888 

± 

0.048 

17.127 

± 

0.022 

0.209 

825 

-2.808 

± 

0.185 

16.981 

± 

0.220 

0.248 

93 

J... 

-3.008 

± 

0.041 

16.748 

± 

0.026 

0.180 

790 

-2.991 

± 

0.144 

16.684 

± 

0.171 

0.195 

93 

H... 

-3.114 

± 

0.040 

16.506 

± 

0.025 

0.173 

781 

-3.149 

± 

0.131 

16.505 

± 

0.156 

0.178 

94 

K... 

-3.203 

± 

0.057 

16.508 

± 

0.040 

0.173 

532 

-3.092 

± 

0.129 

16.364 

± 

0.153 

0.176 

95 

3.6/j-m 

-3.185 

± 

0.034 

16.413 

± 

0.021 

0.146 

789 

-3.259 

± 

0.112 

16.468 

± 

0.134 

0.152 

92 

4.5/im 

-3.166 

± 

0.035 

16.365 

± 

0.023 

0.155 

797 

-3.108 

± 

0.115 

16.276 

± 

0.137 

0.156 

93 

5.8/im 

-3.153 

± 

0.098 

16.324 

± 

0.080 

0.143 

252 

-3.229 

± 

0.109 

16.394 

± 

0.130 

0.145 

90 

8.0/im 

-2.946 

± 

0.909 

16.121 

± 

0.872 

0.146 

42 

-3.254 

± 

0.121 

16.409 

± 

0.143 

0.161 

92 

w... 

-3.320 

± 

0.031 

16.378 

± 

0.019 

0.136 

817 

-3.274 

± 

0.115 

16.329 

± 

0.136 

0.151 

92 


Note. — The subscripts 5 and l, stand for Cepheids with 0.4 < log(P) < 1.0 (i.e. short period Cepheids) and 
log(P) > 1.0 (i.e. long period Cepheids), respectively. ZP and <7 represents the zero point and dispersion of the 
P-L relation, respectively. N is the number of Cepheids used in deriving the P-L relations. 



Fig. 5.— Left Panel: residuals of the fitted P-L relation in Wesenheit function as a function of logarithmic period. Middle Panel: 
distribution of the residuals presented in the left panel. The fitted Gaussian function to the histogram is shown as dashed (red) curve, with 
fitted Gaussian parameters (^, a) given at top of the plot. Right Panel: the qq plot of the residuals. The (red) line represents the case of 
y = X. The theoretical quantiles represent the quantiles if the residuals are drawn from a Gaussian distribution, while the sample quantiles 
are actual quantiles calculated from the residuals. 


summarized in Table 5 for the SMC P-L relations. The 
largest Fg-value for the SMC P-L relations is found to be 
1.17 in the J-band. Therefore, our F-test results strongly 
suggest that the SMC P-L relations are linear from opti¬ 
cal to infrared bands. Our results are further confirmed 
with a non-parametric random walk test (see Koen et al. 
2007; Bhardwaj et al. 2015). A p(i?)-value from the ran¬ 
dom walk test such that p{R) >0.1 would indicate that 
the null hypothesis of cannot be rejected. We found that 
p{R) > 0.3 in all cases with the exception in the i7-band 
P-L relation, which displays a marginal p-value (0.06) 
from the random walk test. We also emphasize that the 
random walk test is non-parametric and does not make 
any assumptions about homoscedasticity or normality of 
residuals. Since the short period Cepheids in 8.0^m-band 
only occupied a small period range, 0.91 < log(P) < 1.0, 
that will cause the determination of the P-L slope to be 
less accurate (—2.946 ± 0.909, see Table 4), we excluded 
the 8.0/im-band P-L relation from our statistical tests. 

5. COMPARISON WITH LMC P-L RELATIONS AND THE 
UNIVERSAL PERIOD-WESENHEIT RELATION 

Since the data sources, algorithms and codes used in 
this work are almost identical to Paper 1, we can compare 
the P-L slopes found in this work for the SMC Cepheids 
to the linear version of the LMC P-L slopes given in 
Paper 1. This provides a critical test for the assumption 
of a universal P-L slopes across different filters. As in 
the previous section, we applied the statistical t-test to 
test the null hypothesis such that P-L slopes from LMC 
and SMC Cepheids be the same in a given band. Table 


TABLE 5 
F-test and 
Permutation Method 
Results for the 
Nonlinearity of P-L 
Relations. 


Band 

Fo 

P 

V... 

0.786 

0.454 

I... 

0.681 

0.504 

J... 

1.165 

0.315 

H... 

0.959 

0.395 

K... 

0.559 

0.580 

3.6/.im 

0.888 

0.414 

4.5/im 

0.526 

0.597 

5.8/im 

0.162 

0.846 

w... 

0.095 

0.908 


6 summarizes our t-test results. We also calculated the 
(absolute) difference of these P-L slopes, |A|, and its as¬ 
sociated error (as), which are also provided in Table 6. 
When comparing the two slopes, it is a common practice 
in the literature to claim the two underlying slopes as 
consistent if A is within, say, ^ 2.5as- In this case, the 
LMC and SMC P-L slopes listed in Table 6 are consistent 
with each other in all bands. Based on the t-test results 
presented in Table 6, the null hypothesis of equivalent 
slopes can be rejected for the P-L slopes in V- and J- 
band only. The i7-band P-L slopes, on the other hand, 
provide evidence of marginal consistency of the slopes. 
Therefore, the slopes of the P-L relations are not univer¬ 
sal in VJ bands, at least for metallicity bracketed by the 
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TABLE 6 

Comparison of the P-L slopes Between LMC and SMC 
Cepheids. 


Band 

Galaxy 

P-L Slope 

|AP 

\T\ 

p-value 

V... 

SMC 

LMC 

-2.660 ±0.040 
-2.769 ±0.023 

0.109 ±0.046 

2.558 

0.011 

I... 

SMC 

LMC 

-2.918 ±0.031 
-2.961 ±0.015 

0.043 ± 0.034 

1.407 

0.160 

J... 

SMC 

LMC 

-3.052 ±0.026 
-3.115 ±0.014 

0.063 ± 0.030 

2.351 

0.019 

H... 

SMC 

LMC 

-3.157 ±0.025 
-3.206 ±0.013 

0.049 ± 0.028 

1.930 

0.054 

K... 

SMC 

LMC 

-3.213 ±0.032 
-3.194 ±0.015 

0.019 ±0.035 

0.616 

0.538 

3.6xtm 

SMC 

LMC 

-3.220 ±0.021 
-3.253 ±0.010 

0.033 ± 0.023 

1.609 

0.108 

4.5^m 

SMC 

LMC 

-3.184 ±0.022 
-3.214 ±0.010 

0.030 ± 0.024 

1.421 

0.155 

5.8^m 

SMC 

LMC 

-3.201 ±0.039 
-3.182 ±0.020 

0.019 ± 0.044 

0.467 

0.641 

S.Oxtm 

SMC 

LMC 

-3.268 ±0.087 
-3.197 ±0.036 

0.071 ± 0.094 

0.805 

0.421 

W... 

SMC 

LMC 

-3.314 ±0.020 
-3.313 ±0.008 

0.001 ± 0.022 

0.055 

0.957 


^ A is the difference of the slopes between LMC and SMC P-L relations. 
The error for A, denoted as < 73 -, is the quadrature sum of the errors in two 
slopes. 


Magellanic Clouds. 

On the other hand, the Wesenheit function, that in¬ 
corporates a color term, provides evidence of a remark¬ 
ably consistent P-L slope between the LMC and SMC 
Cepheids with a value of ~ —3.31. This strongly indi¬ 
cates that the Wesenheit function is universal (for exam¬ 
ple, see discussion in Bono et al. 2010), and hence more 
suitable for the distance scale application than the single 
band P-L relations. Since the P-L slopes for Wesenheit 
function are almost identical, the difference of the P-L 
zero points for LMC and SMC Wesenheit functions di¬ 
rectly translates to the relative distance between LMC 
and SMC. In terms of distance modulus /i, it is found 
to be = 0.483 ± 0.015 mag, which is in good agree¬ 
ment with the value found in Inno et al. (2013, A/x = 
0.48 ± 0.03 mag based on fundamental mode Cepheids 
only) or the preferred value based on Caussian fitting to 
a number of recent measurements (Craczyk et al. 2014, 
Ayx = 0.458 ± 0.068 mag). The recommended distance 
moduli for LMC and SMC are 18.49 ±0.09 mag (de Grijs 
et al. 2014) and 18.96±0.02 mag (de Grijs & Bono 2015), 
respectively, with a difference of A/x = 0.47 ± 0.09 mag. 
Again this value is consistent with our result. By shifting 
the SMC data with 0.483 mag, we combine the Cepheids 
in both Magellanic Clouds and derive the following P-L 
relation for 2578 Cepheids: 

lT = -3.314(±0.009)logP± 15.892(±0.006), a = 0.099, 

which, as expected, is identical to the LMC period- 
Wesenheit relation given in Paper 1. 

6. CONCLUSION 

The main goal of this study was to extend the work 
of Paper I by deriving the multi-band P-L relations for 
SMC fundamental mode Cepheids using the latest com¬ 
pilation of OGLE-111 catalog (Soszyhski et al. 2010). In 
addition to the E/-band mean magnitudes adopted from 
Soszyhski et al. (2010), we also matched the OGLE-111 


SMC Cepheids to 2MASS and SAGE-SMC catalogs and 
derived the mean magnitudes in infrared JiL AT-band and 
the four IRAG bands. Extinction corrections to these 
mean magnitudes were done using the Zaritsky et al. 
(2002) extinction map. These data sources are the same, 
or very similar, to those adopted in Paper 1 for the 
LMC Cepheids. We also applied the almost identical 
algorithms and codes from Paper 1 in this work. The 
difference between this work and Paper 1, on the con¬ 
trary, includes (a) Cepheids with log(P) < 0.4 and some 
longest period Cepheids were excluded in the SMC sam¬ 
ple; and (b) we applied a period cut to the K, 5.8/im- 
and 8.0^m-band P-L relations to avoid the bias due to 
incompleteness at the faint end. 

We then derived the extinction corrected P-L relations 
mVIJHK bands and in the four IRAC bands, as well as 
the extinction free period-Wesenheit function, for SMC 
fundamental mode Cepheids. Following Paper 1, we did 
not apply absolute calibration to our P-L relation be¬ 
cause the readers can adopt their preferred SMC distance 
to calibrate these P-L relations (for example, the latest 
measurement of SMC distance can be found in Graczyk 
et al. 2014). We summarize the main results based on 
the derived SMC P-L relations as follows: 

1. Based on the F-test, the SMC P-L relations are 
found to be linear from optical to infrared bands 
(except the 8.0/xm-band P-L relation at which the 
F-test is not applied) for SMC Cepheids with pe¬ 
riod between log(F) = 0.40 to log(F) = 1.51 (for 
AT-band, the period range is 0.51 < log(F) < 
1.51; for 5.8^m band, the period range is 0.68 < 
log(F) < 1.51). The period-Wesenheit relation is 
also linear in the same period range. 

2. Based on the f-test, the null hypothesis of equiv¬ 
alent slopes for the LMC and SMC P-L relations 
can be rejected in the V J bands. The P-L slopes in 
other bands are consistent between the LMC and 
SMC Cepheids. We also note the remarkable agree¬ 
ment between the SMC/LMC P-L slopes for the 
Period-Wesenheit relations. 

This work is focused on SMC P-L relations, which are 
found to be linear as opposed to the nonlinear LMC P- 
L relations reported in our Paper I (Ngeow et al. 2009) 
and reference therein. The work of Kanbur et al. (2004), 
Kanbur & Ngeow (2006) and Kanbur et al. (2007) has 
provided one possible theoretical scenario by which the 
LMC P-L relation can be nonlinear whilst the SMC P-L 
relation is linear. These relate to metallicity differences 
leading to different mass-luminosity relations obeyed by 
Cepheids in the Magellanic Clouds and a different hy¬ 
drogen ionization front-stellar photosphere interaction in 
terms of periods and phases. The result of identical P- 
L slopes in LMC and SMC period-Wesenheit relations 
is interesting, because the EFband P-L slopes for LMC 
Cepheids are nonlinear but linear in the case of SMC 
Cepheids. In terms of period-color relation, Kanbur & 
Ngeow (2004) found that the mean light period-color re¬ 
lations in LMC and SMC are nonlinear and linear, re¬ 
spectively, but Sandage et al. (2009) suggested the SMC 
period-color relation should be nonlinear. To reconcile 
the linear and identical slopes in LMC and SMC period- 
Wesenheit relations implies that the {V — I) period-color 
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relation must play a substantial role here. We will ad¬ 
dress the issue of period-color relations in a future paper. 
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APPENDIX 

CORRECTION OF T-TEST RESULTS FOR LMC CEPHEIDS BETWEEN OGLE-II AND OGLE-III CATALOGS IN PAPER I 

Due to a mistake made in the code for performing the t-statistical test, the p-value was mis-labelled as expected 
t-value based on a t-distribution given a and v. Therefore, the label “t” in Table 2 of Paper I (Ngeow et al. 2009) 
should be replaced by “p-value”, and hence the LMC P-L slopes based on OGLE-II and OGLE-III catalogs are all 
consistent with each other in all bands. The only exception is the P-L slope for Wesenheit function from Udalski et 
al. (1999b,c, —3.277 ± 0.014 as compared to —3.313 ± 0.008 derived in Paper I). 
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